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Abstract. Regge theory provides an excellent description of small- :r structure- 
function data from Q 2 = up to the highest available values. The large-Q 2 data 
' should also be described by perturbative QCD: the two descriptions must agree in 

. the region where they overlap. However, at present there is a serious lack in our 

understanding of how to apply perturbative QCD at small x. The usual lowest- 
order or next-to-lowest order expansion is not valid, at least not until Q 2 becomes 
^ , much larger than is usually assumed; a resummation is necessary, but as yet we do 

not know how to do this resummation. 

1 Introduction 

Perturbative QCD has become a well-established description of hard pro- 
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cesses. But it is not complete: it must be supplemented with other descrip- 
tions. 

' Regge theory is an example of another description. It was extensively 

C developed some 40 years agoQ and is based on our knowledge of analyticity 

' properties of scattering amplitudes. It relates the high-energy behaviour of 

scattering amplitudes to exchanges of known particles. In order to describe 
the data, it also introduces extra terms that, at least so far, are not related to 
Oh, exchanges of known particles. These extra terms are called pomeron-exchange 

O 1 terms, after the Russian physicist Isaac Pomeranchuk. The notion of the 

D , pomeron is a very old one. We do not know whether it too describes the 

exchange of particles, but if it does there is general agreement that these are 
likely to be glueballs. 

Regge theory gives an excellent description not only of soft hadronic 
, procescs[|][|[j|], but also of the behaviour of the structure function ^(a;, Q 2 ) 

' at small x, all the way from real-photon-induced events (Q 2 = 0) up to very 

deeply inelastic ones (large Q 2 ). The structure function describes an example 
of a semihard process. Another example of a semihard process is the quasi- 
elastic reaction J/ip p. For such semihard processes, the data reveal j| 
that there exists a second pomeron, the "hard" pomeron, in addition to the 
original "soft" pomeron which enters into purely hadronic reactions. The soft 
pomeron is certainly nonperturbative in origin, though it is possible that the 
hard pomeron is associated with the perturbative BFKL equation||. How- 
ever, there are many serious problems^) JtJ with the BFKL equation, and so 
this is not sure. 

Certainly, one would hope that at small x and large Q 2 Regge theory 
can be made to agree with perturbative QCD, in particular with DGLAP 
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evolution|8|. However, there are problems with the DGLAP equation too at 
small x. The equation involves a kernel or splitting function P(z, as{Q 2 ))- If 
one expands this in powers of as, each term is singular at z = 0. HoweverB 
it is rather sure, from general considerations of the known analyticity in Q 2 
of the structure function, that P(z,as{Q 2 )) is not singular at z = 0. The 
singularities in the terms of the expansion are a signal that the expansion is 
illegal near z — 0, and a resummation is needed to remove them. At present, 
we do not know exactly how to do this. It is possible that the presence or 
absence of the singularity is immaterial if Q 2 is large enough fic|| , but for most 
applications in the literature fll|| it is likely to be a real problem. 



2 Regge theory 

Through the optical theorem, the structure function Fi(x,Q 2 ) is the imag- 
inary part of the virtual Compton amplitude Ti(v,t,Q 2 ) evaluated at zero 
momentum transfer, t = 0. Here 2v = 2p.q = Q 2 /x. Regge theory begins|ij 
by considering the crossed-channel process 7*7* — > pp, for which y/i is the 
centre-of-mass energy. It first makes a partial-wave series expansion in this 
channel, in terms of partial- wave amplitudes a^(t, Q 2 ) and Legendre polyno- 
mials Pi(cos8t), where t is the crossed-channel scattering angle. This ex- 
pansion has a definite, but limited, region of convergence in the space of the 
three variables cos9 t ,t,Q 2 or, equivalently, v,t,Q 2 . The partial- wave ampli- 
tude a^(t, Q 2 ) is defined initially for physical values of the angular momentum, 
£ = 0, 1, 2, . . .. By a well-defined procedure, its definition is extended to all 
values of £, both real and complex, and the partial wave series is converted 
to an integral: 



Ti(v,t,Q ) = — / dl : a(£,t,Q ) 



(1) 



The contour C is initially that of figure la, wrapped arround the positive 




(a) (b) 
Fig. 1. Contour C in the complex I plane for the integral (hi) 



real-^ axis on which are located the zeros of the denominator sin n£. But the 
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properties of a(£ , t, Q 2 ) allow one to use Cauchy's theorem to distort it to 
become parallel to the imaginary-i? axis, as in figure lb. It turns out that this 
extends the region of convergence beyond that of the original series, and the 
integral can be continued analytically to where we need it to calculate the 
structure function Fi(x,Q 2 ): 

t = 2v>Q 2 >Q (2) 

As we continue analytically into this region, the singularities of a(£, t, Q 2 ) will 
move around in the complex f-plane, and one of them may try to cross the 
contour C. We must distort the contour again so as to avoid this happening. 
Depending on whether the singularity is a branch point or a pole, we then 
have either figure 2a or 2b. 







5> 



(a) 



(b) 



Fig. 2. Distortion of the contour C caused by (a) a pole, or (b) a branch point, of 
a(£,t, Q 2 ) trying to cross it 



From what is known about the analyticity structure, it is rather sure that 
the positions of the singularities of a(£, t, Q 2 ) in the complex £ plane do not 
depend on Q 2 , only on t. In the case of a pole 

the integration around it yields a contribution 



irl3(Q 2 ,t)P a{t) (-cos6 t ) 
sin na(t) 

to Ti(is, t,Q 2 ). Now when t = the analytic continuation of cos# t is 



(4) 



cos6> t = — — (5) 

iQlTlp 

and when this is large the analytic continuation of the Legendre polynomial 
has the simple behaviour 

-—— £ — (-20089.) (Rc£>- 2 ) (6) 
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So, with (§), we see that the integral ([j]) becomes just a Mellin transform, 
and the "Regge trajectory" a(t) contributes at t = 

MQ 2 K (0) (7) 

to the large-f behaviour of T\(v, 0, Q 2 ). Here, b\(Q 2 ) is a constant multiple of 
/3(Q 2 , 0). In the case of T 2 {v, 0, Q 2 ) its definition includes a kinematic factor 
which reduces the power of v by one unit, so since v = Q 2 jlx this gives 

F 2 (x,Q 2 )~ f(Q 2 )x^ (8) 

with 

e = a(0) - 1 (9) 

Regge theory gives no information about the function f(Q 2 ), other than that 
it is an analytic function with singularities whose locations are known 
The power (1 — a(0)) is independent of Q 2 . 



3 Fit to data 

In the case where the singularity that crosses the contour C is a branch point 
instead of a pole, dragging a branch cut with it as shown in figure 2b, the 
simple power of x in (0) is replaced with something more complicated. One 
knows, from unitarityM, that if there are poles in the complex i plane, there 
must certainly also be branch points. 

On the principle that it is usually the best strategy to try the simplest 
possible assumption first, Donnachie and I tested the hypothesis that at t = 
the contribution from branch points is much weaker than from poles. Not 
everybody agrees with this strategy but we applied this several years 
agoj| to purely hadronic total cross-sections and found that they are all 
described well by a sum of just two powers v ei and v t2 The two powers are 

ei=0.08 ("soft pomeron" exchange) , , 

e 2 = —0.45 (p, uj, f, a exchange) 

More recently Q], we tested the hypothesis that also the contribution to 
the structure function F% at small x from branch points is much weaker than 
from poles. We made a fit of the form 

2 

F 2 (x,Q 2 )~J2fi(Q 2 >- ei (11) 

i=0 

We fixed the values of the powers £i and e 2 to be the same as in ©, and left 
£o as a free parameter to be determined from the small-a: structure-function 
data. Our fitting procedure had three stages. First, we arrived at a provisional 
value for eg by using data only in the region 

x < 0.07 < Q 2 < 10 GeV 2 (12) 
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which gave us 

e w 0.4 (13) 

Next, with this value for eo, we fitted the data to a sum (|ll]) of three powers 
of x for each values of Q 2 for which there exist data, still restricting x to 
less than 0.07 . This gave us the plots shown in figure 3 of the hard and soft 
pomeron coefficient functions fo(Q 2 ) and fi(Q 2 ) as functions of Q 2 . The data 
do not constrain the (/, a)-exchange coefficient function fz{Q 2 ) at all well, 
so we retained the simple form we chose for it in the first stage of our fitting 
procedure. 

We know from gauge invariance that must vanish linearly in Q 2 as 
Q 2 — > at fixed v. We see that the hard-pomeron coefficient function fo(Q 2 ) 
recovers from this rather slowly as Q 2 increases, until it begins at Q 2 w 10 
GeV 2 to rise quite rapidly. On the other hand, the soft-pomeron coefficient 
function fi(Q 2 ) rises rapidly away from Q 2 — 0, until it peaks at between 5 
and 10 GeV 2 . It then falls again. This was a surprise to us: the soft-pomeron 
contribution to the small- X structure function apparently is higher twist. 

The third part of our fitting procedure chooses functions that have the 
general shape of the plots in figure 3, with a number of parameters. It still 
retains the same form for f2{Q 2 ), which matters significantly only for the 
real-photon data at low energy — our fit includes these for J~s > 6 GeV, and 
including the f% term is essential to fit them. So we use ( |ll|) with 

/o(Q 2 )^o(^) 1+e °(l + g)^ (14) 

/ 1 (Q 2 ) = ^ 1 (^) 1+ei (i + V / f)" (15) 

HQ 2 ) =^(crfy 1+£2 (i6) 

The hard-pomeron power is now again a free parameter, together with the 
three coefficients Ai, the mass scales Qf, and Q 2 S . With these 8 parameters, 
we obtain a \ 2 per data point of 1.0 for 595 data points. These data points 
have x < 0.07 and range from Q 2 — to 2000 GeV 2 . Sample plots are shown 
in figure 4. 

I want to make a number of comments on these fits: 



The choice (16) for the analytic form of the coefficient functions is an 
economical set that describes the data well, but there are other choices 
that agree also with the data extracted in figure 3. Using different choices 
results in different values for the hard-pomeron power eo, but they are 
all within 10% of 0.4. 



The choice (16) makes fo(Q 2 ) ~ Q e ° at large Q 2 , which corresponds 
to /3(Q 2 ,0) in (||) becoming constant for large Q 2 . There is no general 
theory that explains this, though it has been predicted ]I3J] from the BFKL 
equation. Parametrisations of fo(Q 2 ) that behave logarithmically at large 
Q 2 can also fit the data satisfactorily. 
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Fig. 3. The hard and soft pomeron coefficient functions fo(Q 2 ) and fi(Q 2 ) ex- 
tracted from the data 



The choice (16) makes fi(Q 2 ) ~ 1/Q at large Q 2 . Fits that make it 
instead behave as 1/Q 2 are also acceptable, though less good. A 1/Q 
behaviour has been predicted p]] from a combination of the BFKL equa- 
tion and infrared renormalons, though this is controversial ]l6| . Fits in 
which fi(Q 2 ) does not go to zero at high Q 2 do not work: it really does 
seem that the soft-pomeron contribution to F2(x, Q 2 ) is higher twist. At, 
say, Q 2 = 5 GeV 2 soft-pomeron exchange dominates in F-z{x, Q 2 ) un- 
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Fig. 4. The Regge fit compared with some of the data, at the largest available Q 
small Q 2 , and Q 2 = 
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til x is less than about 10~ 3 . The consequence of this for conventional 
structure- function fits[ll| needs discussion. 

• Without the contribution from the hard pomeron, the fit to the real- 
photon data agrees wellj^] with the measured HERA points. However, it 
is seen in figure 4 that including the hard pomeron makes the fit pass 
significantly above these. 

• Figure 5 shows a fit to the preliminary ZEUS data for the charm structure 
function. It shows that, to a good approximation, these data may be 
described well by hard-pomeron exchange alone. The fit is actually a 
single-parameter fit; it includes the constraint that at high Q 2 the hard 
pomeron is flavour blind. 

• Figure 6 shows that the data for the semihard process 7p — > J/ipp are 
well described by a mixture of soft and hard pomeron exchange The 
figure is for Q 2 = 0; the HERA data indicate that, as Q 2 increases, the 
hard-pomeron component becomes relatively more important. One might 
guess that the data for 7 p — > pp may be described similarly, though at 
Q 2 = the soft pomeron dominates. 



4 Perturbative evolution 

At high Q 2 , parton distributions evolve with Q 2 according to a perturbative 
evolution equation, called the DGLAP equation. In the singlet channel, the 
equation is written^] in terms of a 2-component vector u whose elements are 
the singlet quark distribution and the gluon distribution: 

f 1 dz 

u(x,t)= —P(z,t)u(x/z,t) (17) 

Jx * 

On the left-hand side, u appears differentiated with respect to 

* = log(Q 2 M 2 ) (18) 

On the right-hand side P is the 2x2 splitting matrix. 

The splitting matrix P may be expanded as power series in as(Q 2 )- It was 
observed long ago|l7| that if one uses just the first term in this expansion, 

p(M) = i^P) pW <19) 

and assumes that each element of u(x,t) has simple power behaviour in x as 
in (||), then the DGLAP equation gives a simple differential equation for the 
coefficient matrix f (Q 2 ) that multiplies the power. Its solution is 

q2\7(«) 



f(Q 2 )=f(^)(log-^) (20) 
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Fig. 5. Hard-pomeron fit to preliminary ZEUS data for F 2 C 




Fig. 6. Two-pomeron fit to data for 7p — ► J/tjjp 
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where "f(N) is the eigenvalue of f3op(N). Here p(N) is the Mellin transform 
of p(JV): 



p(N) = [ dzz N p{z) 
Jo 



(21) 



For e = 0.4 the eigenvalue 7(e) is close to 3. One may fit the data satis- 
factorily by requiring the hard-pomeron coefficient function fo{Q 2 ) to have 
this behaviour at large Q 2 , though not as well as with the power behaviour 
of (|l6|). However, in any case it is not valid to use the approximation ([l9]) at 
low x, even if it is supplemented by adding in higher-order terms in the per- 
turbative expansion. Such an expansion is not legal at small x, even though 
it is widely usedQ. 

The problem is that two of the elements p q c and paa of the matrix p(z) 
have poles[§ at z = 0, so that their Mellin transforms have poles at N = 0. 
For example, 

ttPgg(z) ~ itPgg(N) ~ ^ (22) 

If one includes higher-order terms in the perturbative expansion, all the ele- 
ments of the Mellin transform of the splitting matrix are singular at N = 0. 
The consequence is that, as soon as Q 2 is large enough for the DGLAP 
equation to be applicable, Q 2 > Q\ say, the Mellin transforms of the par- 
ton distributions acquire rather nasty singularities at N = 0, something like 
exp(l/iV). This conflicts with what we believe we know about analyticity 
properties in Q 2 . As I explained in section 2, the Regge amplitude at high 
energy (or small x) is essentially just the Mellin transform: the relationship 
is 

£ = N + 1 (23) 

So a singularity at N = corresponds to one at £ = 1. However, the analyt- 
icity properties of the amplitude whose imaginary part is F2(x, Q 2 ) were well 
studied nearly 40 years ago|lj and there was never any suggestion that it has 
a singularity at £ = 1, let alone a nasty one. If there is no singularity when 
Q 2 is small, it is not possible that one suddenly appears when we continue 
analytically in Q 2 to Q 2 > Q 2 y That is, there can be no singularity of the 
splitting matrix at N = — and probably not, indeed, at any other value of 
N. 

The appearance of the singularity at N = arises because of the pertur- 
bation expansion. Compare, for illustration, the function 



<f>{N,a s ) = N -y/N 2 - as (24) 
Its expansion in powers of as is 

« N > a *) = m + uk + --- (25) 
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Each term in the expansion is singular at N — 0, but manifestly the complete 
function 4>(N, as) is not. The expansion is illegal near N = 0. 

The timelike splitting matrix is very similar [[l8| in form to (j)(N,as), but 
in the spacelike region things are a little more complicated. If one combines 
the DGLAP equation with the BFKL equation, which could well be a valid 
thing to do when x or N is small, one finds that the element Pgg of the 
spacelike splitting matrix is given JTsj] in terms of the Lipatov characteristic 
function g x(u),a s ): 

x(P gg (N,t) 7 a s )=N (26) 

To lowest order in as, 

X(w, a s ) = — pty(l) - " V-(l - «)] (27) 

One may easily verify || that, if one inserts ( |27j ) into (p(f), each term in the 
expansion of P gg (iV, f) in powers of as is singular at N = 0, but the complete 
function P gg (N,t) is finite there. Some authors[20| therefore advocate that 
one should use the solution to the full equation (2q) and avoid expanding it. 

This resummation is obviously sensible, but there are still serious unsolved 
problems. Firstly, it is highly doubtful that it is valid to use the uncorrected 
BFKL equation Q: it assumes that nonperturbative effects cause no complica- 
tions, and does not properly take account of energy conservation. Secondly, 
even if one ignores these difficulties, using the lowest-order approximation 
( p7| ) to x(w,as) is not valid: the next-to-leading-order correction is huge 
There have recently been some interesting attempts to solve this problemf 
but more work remains to be done. 

So at present, although we know that the elements of the splitting matrix 
cannot diverge at N = 0, we cannot calculate them. Any application of the 
DGLAP equation in which they are expanded perturbatively, and in which 
use is made of the fact that the terms in this expansion are large, cannot be 
trusted. 



5 Summary 

• When we use the DGLAP equation at small x it is not valid to use an 
unresummed perturbative expansion of the splitting matrix — at least, 
not until Q 2 is somewhat larger than is normally assumed. At present, 
we do not know how to perform the necessary resummation properly. 

• Regge theory with two pomerons fits the small- a; data for F 2 (x,Q 2 ) ex- 
tremely well up to the highest available values of Q 2 . It strongly suggests 
that at, say, Q 2 = 5 GeV 2 , most of ^(x, Q 2 ) at small x is higher twist. 

• Perturbative QCD and Regge theory are not rival theories. They com- 
plement each other and we have to learn how to make them fit together. 
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